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Strong interactions can amplify quantum effects such that they become important on macroscopic 
scales. Controlling these coherently on a single particle level is essential for the tailored preparation 
of strongly correlated quantum systems and opens up new prospects for quantum technologies. 

Rydberg atoms offer such strong interactions which lead to extreme nonlinearities in laser coupled 
atomic ensembles. As a result, multiple excitation of a Micrometer sized cloud can be blocked while 
the light-matter coupling becomes collectively enhanced. The resulting two-level system, often 
called “superatom”, is a valuable resource for quantum information, providing a collective Qubit. 

Here we report on the preparation of two orders of magnitude scalable superatoms utilizing the 
large interaction strength provided by Rydberg atoms combined with precise control of an ensemble 
of ultracold atoms in an optical lattice. The latter is achieved with sub shot noise precision by 
local manipulation of a two-dimensional Mott insulator. We microscopically confirm the superatom 
picture by in-situ detection of the Rydberg excitations and observe the characteristic square root 
scaling of the optical coupling with the number of atoms. Furthermore, we verify the presence of 
entanglement in the prepared states and demonstrate the coherent manipulation of the superatom. 

Finally, we investigate the breakdown of the superatom picture when two Rydberg excitations are 
present in the system, which leads to dephasing and a loss of coherence. 

PACS numbers: 32.80.Ee, 33.80.Rv, 03.67.Bg, 03.67.-a 


Nonlinearities in light-matter coupling are usually 
weak, leading to a linear growth of the number of optical 
excitations with increasing photon flux. In contrast, the 
most extreme regime of strong nonlinearities is reached 
when an ensemble of many absorbers can host only a 
single excitation, such that one photon already saturates 
the medium. This can be realized with the aid of op¬ 
tical cavities [l] or, in free space, by atomic ensembles 
excited to Rydberg states [2]. For the latter, extremely 
strong dipolar interactions between Rydberg atoms block 
all but a single optical excitation in a volume of several 
Micrometers pu, effectively transforming the N atoms 
within this volume to one collective two-level system. Un¬ 
der uniform illumination this “superatom” features en¬ 
hanced coupling to the light field and the Rydberg ex¬ 
citation is symmetrically shared between the individual 
atoms [7]. The resulting singly excited Dicke state is 
also known as FF-state whose many-body character is 
reflected in multipartite entanglement between its con¬ 
stituents [8]. Superatoms are valuable resources for quan¬ 
tum information. They have been proposed as collective 
Qubits [4] and indeed, strong interactions between them 
were demonstrated recently [9]. These collective Qubits 
have distinct advantages over single atoms that have pre¬ 
viously been entangled using the strong Rydberg inter¬ 
actions mmm- First, the inherent collective enhance¬ 
ment of the atom-light coupling provides a single pho¬ 
ton interface and efficient entanglement transfer between 
atoms and light [HHH]- Second, the information is re¬ 
dundantly stored in the N constituent particles, protect¬ 
ing it against detrimental atom loss mm- Further ap¬ 


plications reach from advanced Qubit encoding schemes 
in multi-level atoms da m to efficient single photon 
sources duns and single photon subtraction [21]. Co¬ 
herent manipulation of superatoms is at the heart of these 
proposals and amounts to controlling the strong and spa¬ 
tially dependent interactions of Rydberg atoms, which, in 
larger samples, lead to dephasing and prohibit the clear 
observation of Rabi oscillations [22H29] . However, for 
small systems of up to 16 atoms Rabi oscillations have 
been directly observed [6l l30lf32] , while for larger systems 
indirect detection methods were required m- In con¬ 
trast, we studied isolated superatoms, locally controlling 
the shape on the single atom level with sub-Micrometer 
precision and the atom number fluctuations better than 
shot noise. We detected the Rydberg excitations in-situ 
with single atom sensitivity and coherently manipulated 
collective systems with scalable size between one and 185 
individual atoms. In this work we report on these mea¬ 
surements and confirm the predicted y/N enhancement 
of the Rabi coupling over two orders of magnitude. Ad¬ 
ditionally, we detect entanglement between the compo¬ 
nents of the superatom and analyze the local distribution 
of the Rydberg excitation within the ensemble. Finally, 
we show that multiple excitations indeed lead to dephas¬ 
ing which, however, is reduced in the post-selected single 
excitation subspace. 

Our superatoms were formed out of an ensemble of ul¬ 
tracold atoms held in a two-dimensional optical lattice 
with unity occupation per lattice site [29] . This system 
was then approximately uniformly coupled to a Rydberg 
state with coupling strength Q. The atoms occupied the 
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FIG. 1. Superatom preparation, (a) Illustration of the symmetric ground and singly excited state (VF-state). Left: iV-atom 
collective Bloch sphere with excitation numbers and coupled states (south pole and IF-state, represented by the red plane) 
indicated. The small pictograms above and below the sphere depict the lattice system with atoms in the ground (red) and 
Rydberg (blue) state. The Husimi distribution of these states and their enhanced coupling Qn is shown in the center. This 
accessible state space defines a superatom represented by the standard Bloch sphere on the right, (b) Atom number histograms 
of the initially prepared samples (blue bars) with Gaussian fits (solid green line). The numbers give the mean and standard 
deviation (s.d.) for each data set. Measured and reconstructed occupation of lattice for exemplary initial states is depicted 
above the respective histograms, c.f. schematic pictograms in (a). The Poissonian distribution with the same mean atom 
number is shown as a reference (green dashed line), (c) Averaged initial ground state atom distributions for the respective 
histograms above. The size of blockade radius Rb is shown by the blue bar. (d) Rabi oscillation data (blue points) and sinusoidal 
fits with exponentially decaying contrast (solid gray line) for iV=7.7(2.2) and iV=131(7). The red line shows the same fit on 
an axis scaled to the number of ground state atoms N (right axis). All errorbars are s.e.m. 


Rydberg state only for a few Microseconds, such that 
their motion in the optical lattice, typically on a Millisec¬ 
ond time scale, could be safely neglected. The excited 
state was chosen such that for most of the experiments 
presented here the system size was much smaller than 
the dipole blockade radius R 5 . This dipole blockade origi¬ 
nates from the van-der-Waals interaction which causes an 
energy shift hA v aw = between Rydberg atoms sep¬ 


arated by the distance R [4 . The extraordinary strong 
interaction tunes the excitation laser out of resonance at 
the blockade distance Rb = (^) such that the system 
is restricted to a single Rydberg excitation. This single 
excitation is symmetrically shared among all N atoms 
if both coupling and interaction are effectively uniform. 
Hence, the system dynamics is confined to the symmetric 
subspace of zero (n e = 0 ) and one (n e = 1 ) excitations, 
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FIG. 2. Collective Rabi oscillations, (a) Collective Rabi oscillation data of the mean Rydberg atom number N e (blue points) 
for different numbers of ground state atoms N = 185(8), 84(6), 42(4), 20(3), 0.74(0.60) (top to bottom) with exponentially 
decaying sinusoidal fits (gray). All errorbars are s.e.m.. (b) Density of detected Rydberg atoms for the datasets in (a) with 
normalized vertically or horizontally averaged density (blue solid line) compared to the initial state atom distributions (red 
solid line). (c) Histograms of the Rydberg excitation number integrated over the total oscillation (blue bars) and at the position 
of the first maximum (orange bars, position in (a) marked by orange solid line). 


whose basis are the Fock states |0) = |gi,..., g^) and the 
entangled TT-state |1) = ^ 7 = • • • , r i, ■ ■ ■ ,9 n), 

where gi and r* label the i-th atom in the ground or 
Rydberg state. Then, the Hamiltonian can be writ¬ 
ten in the simple form H = hy/N Q / 2(|0)(1| + |1)(0|), 
where the symmetry induced collectively enhanced cou¬ 
pling y/~NQ appears explicitly. The collective Bloch 
sphere (Fig. 00 offers a convenient way of representing 
states within the symmetric subspace via their Husimi 
quasiprobability distribution [331 EH- The state |0) lies 
at the south pole of the sphere, while |1) corresponds to 


a ring-like structure. These many-body states form the 
basis for the superatom. 

For the preparation of the superatoms our experi¬ 
ment started with a two dimensional degenerate gas of 
rubidium-87 in the |F = 2 ,rriF = 2) hyperfine state, con¬ 
fined in a single antinode of a vertical (z- axis) optical 
lattice at 80 E r . Here, E r is the lattice recoil energy of 
our square lattice with periodicity a\ at = 532 nm. We 
prepared a unity filling Mott insulator of ~ 200 — 500 
atoms by adiabatically switching on two orthogonal lat¬ 
tices in the x — y plane to 40 E r . We then used our lo- 
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Ground state atom number N 

cal addressing technique to precisely control the size and 
shape of the atomic sample to a square with diagonal 
length D containing between one and 185 atoms [35, 36 . 
This ensured that the edges and the total atom number 
of the atomic samples were well defined, allowing us to 
measure total atom number fluctuations up to 4 dB below 
shot noise (Fig. []J),c). The atoms were then coupled to 
the 68S'i/ 2 , | mj = —1/2) Rydberg state via a two-photon 
scheme (red laser with wavelength 780 nm and blue laser 
with wavelength 480 nm) [29 . The excitation beams were 
counterpropagating perpendicular to the atomic plane 
(^-direction) with waists wo = 44(2)/im for the red and 
wo = 17(5)/im for the blue beam. The van-der-Waals 
coefficient of the 6SS state is Cq = h • 630 GHz /i m 6 , re¬ 
sulting in blockade radius of R 5 = 11.7(1) fim for the 
single particle Rabi frequency of £2 = 2tt • 240(11) kHz. 

We detected the Rydberg atoms in-situ using an ef¬ 
ficient (> 99.9%) push-out of the ground state atoms 
lasting 8/is followed by a stimulated depumping of the 
Rydberg atoms back to the ground state. The remain¬ 
ing atoms were then imaged using in-situ fluorescence 
detection with a position resolution of approximately ±1 
lattice site [29] . From our data we inferred an overall effi¬ 
ciency of rj = 0.67(5) for the spatially resolved detection 
of a single Rydberg atom. The spatial control over the 
sample allowed for microscopic control of the superatom 
size. As long as we ensured D (up to and includ¬ 

ing 131 atoms) we observed coherent enhanced Rabi os¬ 
cillations between the zero and one excitation subspaces. 
Here, the ratio of the amplitude of the Rabi oscillations 
to the total atom number scales as 1/N as opposed to 
being constant for independently oscillating particles. In 
Fig.[p we illustrate this scaling for two exemplary cases 
of 8 and 131 atoms. 

In order to characterize the prepared superatoms mi- 


FIG. 3. Scaling of the Rabi frequency and entanglement, (a) 
Extracted values of £2n (blue points) versus average initial 
atom number N for the data shown in Figs. 0 and [2] with 
a power law fit (green line). The inset shows the exponen¬ 
tial decay time of the Rabi oscillations (blue points). The 
expected decay based on the reference sample atom num¬ 
ber fluctuations (dark green shading) and, additionally, tak¬ 
ing into account noise in the pulse area (light green shading) 
are shown for comparison, (c) Extracted visibility V of the 
collective Rabi oscillation versus atom number N after one, 
three and five half cycles of oscillation (red data points with 
increasing lightness, shifted slightly horizontally for better 
visibility). The blue shaded area includes all classical states 
with fully separable single particle density matrices. The gray 
points show the visibility after one half cycle corrected for the 
measured detection efficiency. All errorbars are la statistical 
uncertainty from the fits. 


croscopically, we measured the spatial distribution of the 
observed Rydberg atom and the excitation statistics dur¬ 
ing the Rabi oscillation (Fig. [2|. For different sample 
sizes between one and 185 atoms we drove Rabi oscilla¬ 
tions by illuminating the sample with the coupling lasers 
for varying duration T. For each T we repeated the ex¬ 
periment 25 — 30 times and extracted the mean Rydberg 
number N e (Fig. [ 2 ^ 1 ). The dramatic acceleration of the 
Rabi oscillation with N is clearly visible in the data. 
Additionally, we compare the spatial distribution of the 
Rydberg atoms (integrated over all T) to the initial dis¬ 
tribution of ground state atoms. Within statistical uncer¬ 
tainty we find a flat distribution consistent with the uni¬ 
form coupling assumption (Fig. J 2 Jd). We experimentally 
confirm the picture of a fully dipole blockaded sample by 
extracting the histogram of the Rydberg excitation num¬ 
bers n e both integrated over the whole observation time 
T and as well at the 7 r-pulse time T n . For sample sizes 
up to 131 atoms, the probability of measuring doubly ex¬ 
cited states with two detected Rydberg atoms was below 
1%. We obtained typically 1 — 4 images with two exci¬ 
tations per 500 — 800 shots. For the largest sample used 
in our experiments, the blockade starts to break down 
and the probability increased to 4. 8 (1.0)% (27 events per 
564 shots). None of the data shown here was corrected 
for the detection efficiency and the measured excitation 
number after T n is consistent with N e = 1 when taking 
it into account. 

One striking signature of the superatom is its 
symmetry-enhanced coupling to the radiation field. We 
extracted the oscillation frequency £2^, the decay time r 
and a global offset A from the data shown in Fig. [ljl and 
[2^1 via a fit to N e = 77 • (A — e~ t l T cos (£2 N ' t ) / 2 ) . Indeed, 
we confirm the expected scaling £2n oc \/ r N over two or¬ 
ders of magnitude (Fig. 0 ). A power law fit of the form 
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i? 7 V = ft * N a yields an exponent of a = 0.48(10). Devia¬ 
tions towards higher Rabi frequencies for small N might 
be due to a residual detuning of the coupling lasers that 
we calibrated via spectroscopy on a dilute atomic cloud 
with an uncertainty of ±200 kHz. A systematic lower 
Rabi frequency at large N can be caused by a residual 
inhomogeneity of the laser coupling (up to 10 %) due to 
its Gaussian intensity profile. Also, the observed onset 
of a breakdown of the blockade for our largest prepared 
samples results in a deviation from the v^V-scaling. The 
latter effect is additionally visible in the extracted steady 
state mean Rydberg atom number r]A. For all but the 
largest sample size we find A = 0.51(2), consistent with 
the expected value of 0.5 . For N = 185 it increases 
significantly to A = 0.65(12). To answer the question 
whether the collective speedup can be exploited for quan¬ 
tum operations or whether decoherence effects dominate, 
we analyzed the quality factor of the Rabi oscillations, 
which is given by the product of the decay time r of the 
measured oscillations and the Rabi frequency f2jy. In¬ 
deed, we find a peaking quality factor for N = 131 due 
to the increasing oscillation frequency but constant decay 
time r « 1 /is in the fully blockaded regime (Fig. in¬ 
set). Among the limiting factors for the coherence time 
are residual atom number and coupling power fluctua¬ 
tions (8(2)%). However, these alone cannot explain the 
observed decay (Fig. If inset). For small atom numbers 
additional decoherence might be due to phase noise and 
slow frequency drifts of the lasers, while the inhomogene¬ 
ity in the Rabi coupling becomes significant at larger N. 

Next to the collective enhancement of the optical cou¬ 
pling, the structure of the excited state itself bears the 
marks of the strong particle correlations. The unam¬ 
biguous proof that the excited state of the superatom is 
indeed the N particle entangled VF-state would require 
full state tomography 133, which is not feasible in our 
setup. However, we will show that the experimental ob¬ 
servations are incompatible with the expectations for a 
fully separable state. To this end, we employ the en¬ 
tanglement witness developed in Ref. [ 8 ] which compares 
the overlap of the VF-state and the prepared state to the 
maximally possible overlap with any separable state. In 
order to extract the VF-state overlap from the measured 
oscillation we model the dynamics by a sum of density 
matrices, p meas = ^(t)p c (t) ± P(t)pm ? where the proba¬ 
bilities a(t) and fi(t) interpolate between the coherently 
oscillating part described by p c {t) and a mixed density 
matrix p m showing no dynamics. Under the experimen¬ 
tally verified assumption that maximally one excitation 
was present in the system, p c (t) describes the subspace 
spanned by the ground and the VF-state, as population 
of all other singly excited states leads to dephasing of 
the oscillations of the uniformly coupled system. Hence, 
a(t) can be identified with the overlap of the prepared 
state with the VF-state and it is directly measured by 
the visibility V(t) of the Rabi oscillations. In Fig. we 



FIG. 4. Superatom Ramsey spectroscopy. Evolution of the 
mean Rydberg atom number N e (blue points) versus hold 
time Tr between the two 7r/2-pulses (schematic in the inset) 
and sinusoidal fit with Gaussian decay (gray line) for an initial 
sample size of N — 38(3) atoms. All errorbars are s.e.m. 

show that the extracted V at the oscillation maximum 
after the first half Rabi cycle is consistently above the 
threshold for entanglement, even without correcting for 
the detection efficiency 77 , while it decays into the clas¬ 
sically allowed region for longer times. If the detection 
efficiency 77 is taken into account, we obtain a lower bound 
for the average VF-state preparation fidelity of 0.77(20) 
at the first maximum. 

To investigate the coherence of the collective Qubit 
further, we used a N = 38(3) atom ensemble for Ram¬ 
sey interferometry [9]. First, we prepared a coherent 
superposition of |0) and |1) by a 7 r/ 2 -pulse of length 
T = 7 r/( 2 fiy/N). After a variable hold time Tr, a sec¬ 
ond 7 r/ 2 -pulse was applied and the mean Rydberg atom 
number N e was measured (Fig. [4]). Due to the ac-Stark 
shift created by the red 780 nm laser during the excitation 
pulse, the calibrated transition frequency differs from the 
bare ground to Rydberg state transition frequency, which 
defines the reference for the Ramsey interferometer. 

Therefore, the observed phase accumulation rate is 
given by this ac-Stark shift and agrees well with an in¬ 
dependent calibration of the latter via microwave spec¬ 
troscopy. The extracted decay time of the Ramsey fringe 
of tr = 2.2(4)/is exceeds the damping r of the Rabi oscil¬ 
lations, indicating that there laser power fluctuations and 
coupling inhomogeneities are likely the main source of de¬ 
coherence. We also performed Ramsey interferometry for 
different N and find within experimental uncertainty tr 
approximately independent of N. 

When increasing the sample to a size where the maxi¬ 
mal distance between two atoms D approaches the block¬ 
ade radius, the isolated superatom-picture is expected to 
break down [38] . The gap to doubly excited states with 
two atoms pinned to the diagonal corners of the pre¬ 
pared square shaped density distribution is smallest, such 
that these are the first doubly excited states populated. 
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FIG. 5. Breakdown of the block¬ 
ade. Measurement of the collective 
Rabi oscillation in sample with N = 
185(8) atoms, (a) The contribution 
of states with n e — 1 (blue points) 
and n e = 2 (red points, shifted 
slightly horizontally for better vis¬ 
ibility) agree with the theoretical 
calculation (green and orange solid 
lines, scaled by extracted detection 
efficiency 77 ). For comparison we 
show the fit (gray line) to the mean 
(c.f. Fig. [2b, which shows twice 
faster dephasing than the n e = 1 
subspace alone. The inset shows 
two exemplary single shots with 
n e — 2 (field of view 14 x 14 /im, 
indicated by gray bar) and the two- 
dimensional pair-correlation func¬ 
tion Cij of all n e — 2 events. Color 
scale: red (anti-correlation) to blue 
(correlation). Data binned 4x4 
sites. All errorbars are s.e.m. 


We discussed already several indications of this block¬ 
ade breakdown for the N = 185 atom sample. Here, the 
maximum separation of two atoms was D = 9.8(7) \± m, 
close to the blockade radius R 5 = 11.7 fi m. In Fig. [5] we 
study the effects of the doubly excited states on the decay 
of the Rabi oscillation for this setting by post-selecting 
the acquired data to single and double Rydberg events. 
The decay time of the singly excited component is two 
times larger compared to the full sample and agrees with 
the prediction of a numerical calculation of the system 
in a reduced Hilbert space [29]. This shows that the ob¬ 
served decay is significantly influenced by the dephasing 
due to double excitation. We observe a slow increase 
of the doubly excited fraction (also in agreement with 
theory) that is consistent with the picture of two inter¬ 
acting excitations [39] . Their interaction energy corre¬ 
sponds to the energy shift Z\ v aw at the distance D, re¬ 
sulting in a detuned optical coupling. At the same time 
the collective enhancement of the coupling to this state 
is only y/2 reflecting the two possible orientations along 
the diagonals. The resulting time scale n/y/2 i? 2 + zij dw 
matches roughly the observed slow rise of the doubly ex¬ 
cited states, however, explaining their probability quan¬ 
titatively requires a more complex model including also 
atoms close to the corners of the square. Spatial correla¬ 
tion measurements confirm the localization of the dou¬ 
bly excited events at the diagonal corners (Fig. [ 5 ] in¬ 
set). Low statistics requires here the binning of 4 x 4 
lattice sites for the evaluation of the two point correla¬ 
tion Cij = {{P(x,y)P(x+i,y+j)) — (P(x,y)) (P(x+i,y+j)))x,y 
Here, P(ij) is the probability of finding a Rydberg exci¬ 
tation in bin (i,j) and (-) x ,y and (•) denote the spatial 
and ensemble averages. 

In conclusion, we demonstrated coherent control 


and two-orders of magnitude scalability of Rydberg- 
superatoms. Using in-situ microscopical detection of the 
Rydberg atoms we confirmed the superatom picture and 
proved presence of entanglement in the involved singly 
excited many-body states. We also demonstrated that 
the collectively enhanced coupling can be harnessed to 
increase the fidelity of collective Qubit rotations under 
realistic experimental conditions. The experiments con¬ 
firmed that coupling to many-body states with larger 
Rydberg occupation leads to interaction induced dephas¬ 
ing, strongly supporting the coherent description of our 
previous experiment on short time scales [29, 36!. Our 
results pave the way towards the controlled step-wise 
preparation of higher Dicke states [4], which have been 
proposed for metrology at the Heisenberg limit m, and 
they promise to shed light on macroscopic entangled 
quantum systems [4TJ. 
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